home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlasd8.z / dlasd8
Encoding:
Text File  |  2002-10-03  |  5.2 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSDDDD8888((((3333SSSS))))                                                          DDDDLLLLAAAASSSSDDDD8888((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASD8 - find the square roots of the roots of the secular equation,
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE DLASD8( ICOMPQ, K, D, Z, VF, VL, DIFL, DIFR, LDDIFR, DSIGMA,
  13.                         WORK, INFO )
  14.  
  15.          INTEGER        ICOMPQ, INFO, K, LDDIFR
  16.  
  17.          DOUBLE         PRECISION D( * ), DIFL( * ), DIFR( LDDIFR, * ),
  18.                         DSIGMA( * ), VF( * ), VL( * ), WORK( * ), Z( * )
  19.  
  20. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  21.      These routines are part of the SCSL Scientific Library and can be loaded
  22.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  23.      directs the linker to use the multi-processor version of the library.
  24.  
  25.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  26.      4 bytes (32 bits). Another version of SCSL is available in which integers
  27.      are 8 bytes (64 bits).  This version allows the user access to larger
  28.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  29.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  30.      only one of the two versions; 4-byte integer and 8-byte integer library
  31.      calls cannot be mixed.
  32.  
  33. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  34.      DLASD8 finds the square roots of the roots of the secular equation, as
  35.      defined by the values in DSIGMA and Z. It makes the appropriate calls to
  36.      DLASD4, and stores, for each  element in D, the distance to its two
  37.      nearest poles (elements in DSIGMA). It also updates the arrays VF and VL,
  38.      the first and last components of all the right singular vectors of the
  39.      original bidiagonal matrix.
  40.  
  41.      DLASD8 is called from DLASD6.
  42.  
  43.  
  44. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  45.      ICOMPQ  (input) INTEGER
  46.              Specifies whether singular vectors are to be computed in factored
  47.              form in the calling routine:
  48.              = 0: Compute singular values only.
  49.              = 1: Compute singular vectors in factored form as well.
  50.  
  51.      K       (input) INTEGER
  52.              The number of terms in the rational function to be solved by
  53.              DLASD4.  K >= 1.
  54.  
  55.      D       (output) DOUBLE PRECISION array, dimension ( K )
  56.              On output, D contains the updated singular values.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAASSSSDDDD8888((((3333SSSS))))                                                          DDDDLLLLAAAASSSSDDDD8888((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      Z       (input) DOUBLE PRECISION array, dimension ( K )
  75.              The first K elements of this array contain the components of the
  76.              deflation-adjusted updating row vector.
  77.  
  78.      VF      (input/output) DOUBLE PRECISION array, dimension ( K )
  79.              On entry, VF contains  information passed through DBEDE8.  On
  80.              exit, VF contains the first K components of the first components
  81.              of all right singular vectors of the bidiagonal matrix.
  82.  
  83.      VL      (input/output) DOUBLE PRECISION array, dimension ( K )
  84.              On entry, VL contains  information passed through DBEDE8.  On
  85.              exit, VL contains the first K components of the last components
  86.              of all right singular vectors of the bidiagonal matrix.
  87.  
  88.      DIFL    (output) DOUBLE PRECISION array, dimension ( K )
  89.              On exit, DIFL(I) = D(I) - DSIGMA(I).
  90.  
  91.      DIFR    (output) DOUBLE PRECISION array,
  92.              dimension ( LDDIFR, 2 ) if ICOMPQ = 1 and dimension ( K ) if
  93.              ICOMPQ = 0.  On exit, DIFR(I,1) = D(I) - DSIGMA(I+1), DIFR(K,1)
  94.              is not defined and will not be referenced.
  95.  
  96.              If ICOMPQ = 1, DIFR(1:K,2) is an array containing the normalizing
  97.              factors for the right singular vector matrix.
  98.  
  99.      LDDIFR  (input) INTEGER
  100.              The leading dimension of DIFR, must be at least K.
  101.  
  102.      DSIGMA  (input) DOUBLE PRECISION array, dimension ( K )
  103.              The first K elements of this array contain the old roots of the
  104.              deflated updating problem.  These are the poles of the secular
  105.              equation.
  106.  
  107.      WORK    (workspace) DOUBLE PRECISION array, dimension at least 3 * K
  108.  
  109.      INFO    (output) INTEGER
  110.              = 0:  successful exit.
  111.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  112.              > 0:  if INFO = 1, an singular value did not converge
  113.  
  114. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  115.      Based on contributions by
  116.         Ming Gu and Huan Ren, Computer Science Division, University of
  117.         California at Berkeley, USA
  118.  
  119.  
  120. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  121.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  122.  
  123.      This man page is available only online.
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.